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The Hubbard model on the fee lattice is studied in the limit of infinite spatial dimensions. At suffi- 
ciently strong interaction finite temperature Quantum Monte Carlo calculations yield a second order 
phase transition to a highly polarized ferromagnetic state. Ferromagnetism is stable over a wide 
hole doping regime but not for electron doping. A possible non-Fermi liquid behavior is discussed. 
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More than 30 years ago, the Hubbard model was intro- 
duced to describe correlation effects in transition metals, 
in particular the band ferromagnetism of Fe, Co and Ni. 
However, it appeared to be a generic model rather for 
anii-ferromagnetism than for ferromagnetism. In fact, 
our knowledge on the possibility of ferromagnetism in 
the Hubbard model is still very limited. Only in a few 
special cases was the Hubbard model proven to have fer- 
romagnetic order. The first rigorous result showing a 
fully polarized ferromagneticground state, the theorem 
by Nagaoka and Thouless [[y||] , is only valid for a single 
hole added to the half filled band in the limit of infinite 
on-site repulsion. Its thermodynamic relevance is un- 
clear until today. In the opposite limit of low electronic 
density a fully polarized ferromagnetic ground state has 
been found recently by Muller-Hartmann || in the Hub- 
bard model on a zigzag chain in the continuum limit. 
Further, the existence of a ground state with net polar- 
ization has been proven for the half filled band case on 
bipartite lattices with asymmetry in the number of sites 
per sublattice Q, and in 'flat-band' systems which have 
a ferromagnetic ground state for any interaction strength 
U larger than zero 

Since ferromagnetism in the Hubbard model has been 
found only in some particular cases one might suspect 
that the model is just too oversimplified to explain band 
ferromagnetism even in a qualitative way. One missing 
feature relevant in real materials is band degeneracy. A 
mechanism similar to Hund's rule in individual atoms 
could lead to a ferromagnetic ordering of the electronic 
spins. In addition, all non-local interactions are neglected 
in the Hubbard model. The question if nearest-neighbor 
interactions, in particular the Coulomb exchange and a 
pair hopping term, are able to induce ferromagnetism has 
been adressed by Hirsch and Tang 0,^]. Within mean- 
field approximations and Monte Carlo simulations in one 
dimensions at half filling they found fully and partially 
polarized ground states. Recently, Strack and Vollhardt 
derived rigorous criteria for the stability of saturated 
ferromagnetism in an extendend single band Hubbard 
model including all nearest neighbor interactions at half 
filling djl^]. Again, the exchange interaction plays an 
important role in stabilizing ferromagnetism. 

However, since the results are valid only for the insu- 



lating half filled band case, the question about the neces- 
sary microscopic ingredients for metallic ferromagnetism 
remains open. 

In this paper we are concerned with the 'pure' Hub- 
bard model, i.e. a single electronic band with only on-site 
interaction. The question we want to answer is: 

Is the pure Hubbard model capable of describing 
metallic ferromagnetism? 

Within the pure Hubbard model we still have freedom 
in the choice of the lattice geometry. We choose the fee 
lattice for the following reasons: i) It is not bipartite. 
Thus we expect the competing antiferromagnetic order 
to be suppressed or at least weakend by frustration, ii) 
The fee lattice has been shown to be a good 'environ- 
ment' for ferromagnetism by variational calculations of 
the stability of the Nagaoka state against doping [|ll]-jl3j, 
by exact diagonalization studies |l3|,[l4| , and also by na- 
ture since Ni with fee structure is fully polarized whereas 
Fe (bec-structure) is not, as pointed out in ref. Since 
the Hubbard model cannot be solved in finite dimen- 
sions, d > 1, and Monte Carlo simulations would allow 
only very small fee lattices we study it in the limit of 
infinite dimensions, d = oo Jl5{ ]. Here, the system is 
reduced to a dynamical single site problem |16| equiva- 
lent to an Anderson impurity model complemented by a 
self-consistency condition fl7|| . Still, it cannot be solved 
analytically, and to avoid further approximations we em- 

ta finite temperature Quantum Monte Carlo method 
This approach has proven to be a powerful and re- 
liable tool for the study of correlated fermi systems both 
for thermodynamical as well as dynamical properties [fl9[ . 
Spectral properties are obtained by analytic continuation 
of the imaginary time data using the Maximum Entropy 
method Q. 

A non-trivial generalization of the fee lattice in d di- 
mensions pM is the set of all points with integer cu- 
bic coordinates summing up to an even integer. It is a 
non-bipartite Bravais lattice for any dimension d > 2. 
Each point has Z = 2d(d — 1) nearest neighbors de- 
fined by all vectors x which can be written in the form 
x = ±ej ± ej, with two different cubic unit vectors ei 
and ej — d). With the proper scaling of the 
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hopping term, t = 1/ \Z |l5j] , the density of states of the 
non-interacting electrons can be calculated in d — oo |2l| l 
to be: 



N (E) = e-^+^/VvQl+V^B) 



(1) 



where we choose contrary to usual convention a posi- 
tive hopping integral which gives a strong square-root 
divergency at the lower band edge, — 1/V2, and no up- 
per band edge pjfl . The lattice structure enters in the 
d = oo-approach only via the density of states. The non- 
bipartite structure of the fee lattice leads to the strong 
asymmetry in Nq(E). 

The variational treatment of the Nagaoka state on the 
fee lattice [Tl|-[il| provides upper limits for the critical 
doping 5 C and lower limits for the critical interaction 
U c where saturated ferromagnetism becomes unstable. 
While for the simple cubic lattice in d — oo the stabil- 
ity regime shrinks to the point U c — oo, 6 C = |p3| , for 
the fee-density of states (1) the Nagaoka state is stable 
above a critical line U c (n) with U c (0) = and U c (l) = oo 
p4j . The Nagaoka state is unstable in the case of electron 
doping (n > 1). Antiferromagnetism is not expected on 
the fee-lattice in high dimensions because the difference 
of the numbers of not frustrated bonds and frustrated 
bonds is only of the order of d resulting in an effective 
field of the order of t 2 d oc 1/d || . 

To detect a ferromagnetic instability we calculate the 
corresponding susceptibility, \f, from the two-particle 
correlation functions |25|. We choose a relatively strong 
interaction of U = 4 ]26| and an electronic density of 
n = 0.58, where the ferromagnetic response appeared 
to be strong. As a function of temperature, we find 
Xf to obey a Curie- Weiss law and extrapolate the tran- 
sition temperature linearly from the zero of Xf 1 ( see 
Fig. 1). The numerical result is a Curie temperature 
of T c = 0.051(2) (27). We also calculate the magnetiza- 
tion, M, within the ferromagnetic phase, i.e. below T c . 
M grows strongly at lower temperatures reaching more 
than 80% of the fully polarized value (M max = n = 0.58) 
at the lowest temperature, T = 0.0357. A fit through 
the three data points with a Brillouin function (Fig. 1) 
would give a full polarization at T — and the same 
critical temperature of T c = 0.05 as above. Translated 
to the three dimensional fee lattice with Z = 12 nearest 
neighbors and a bandwidth of W — 16i the critical tem- 
perature becomes T C (3D) ss 0.011VF. Thus, despite the 
oversimplifications of the single band Hubbard model, 
the resulting Curie temperature has a realistic order of 
magnitude of 500-800K for typical values of W around 5 
eV. 

In order to answer the question if the system is still 
metallic we calculate the single particle spectrum. In 
the ferromagnetic phase we can distinguish two differ- 
ent spectra for spins parallel and anti-parallel to the net 
magnetization (Fig. 2). The latter (minority spin) spec- 
trum (dotted line) is shifted to higher frequencies and 
has a pronounced upper band around ui — ^ ~ U = 4 



with a developing (pseudo) gap. The majority spin spec- 
trum, however, is only slightly affected by the interaction. 
Here, the weight of the upper band is very small since 
due to the Pauli principle the majority spin-electrons are 
unlikely to hop over preoccupied sites. Apparently the 
system is metallic since both spectra have a finite value 
at the Fermi level (ui — /1). This also holds in the para- 
magnetic phase (not shown). 
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FIG. 1. Magnetization, M (diamonds), and inverse ferro- 
magnetic susceptibility, Xf- 1 (circles; values multiplied by a 
factor of 4 to use the same scale) for U — 4 and n = 0.58. 
Error-bars are of the size of the symbols or smaller. Dotted 
line is a linear fit to xi? 1 , dashed line a fit with a Brillouin 
function to M. 
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FIG. 2. Single particle spectrum, N(ll>), for (7 = 4, 
n = 0.58, and T = 0.04. Full (dotted) line: majority (mi- 
nority) spin spectrum. 

Seemingly our finding violates the general theorem on 
local Fermi liquids (LFL) [^8| that prohibits ferromag- 
netic instabilities in a single band Fermi liquid with local 
self energy. How can this contradiction be resolved? In 
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the present situation the self energy E„ = E(iw n ) is local 
by construction, therefore we investigate its imaginary 
part to determine if the system is a Fermi-liquid. In that 
case the imaginary part of £„ vanishes at least linearly 
with Matsubara frequency \u> n \ — > 0. Although our finite 
temperature approach does not allow to reach the limit 
w n — * 0, the results close above and below the transition 
(Fig. 3) strongly suggest that Im £ n approaches a finite 
value at low frequencies or at least vanishes much slower 
than uj n . Together with results for the single particle 
spectrum we conclude that the system is likely to be a 
metal with non-Fermi liquid properties. The existence of 
a metallic non-Fermi liquid state in high dimensions is 
quite surprising but has been observed before in an ex- 
tended Hubbard model [^| and also in an alloy analogy 
approach to the single band Hubbard model over a wide 
range of doping and interaction strength ffl] . It can also 
not be ruled out by existing 'essentially exact' treatments 
of the Hubbard model on infinite dimensional bipartite 
lattices |l9[ ]. Further investigations on the nature of the 
paramagnetic phase are required to settle this question. 
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FIG. 3. Imaginary part of the self energy for U — 4 and 
7i = 0.58 in the paramagnetic phase (T = 1/13,1/15,1/20, 
dashed line; curves are on top of each other) and in the ferro- 
magnetic phase (T = 1/25, solid lines; smaller absolute values 
correspond to the majority spin) . Statistical errors are within 
the width of the lines. 



also checked for phase separation by calculating the elec- 
tronic compressibility but never observed any instability. 

The phase diagram in Fig. 4 is qualitatively similar to 
the results of a 'self consistent moment approach' to the 
the Hubbard model on the three dimensional fcc-lattice 
fjlf . In particular the optimal electron density is found 
to be n w 0.65 which is consistent with our results. 
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FIG. 4. Density dependence of the Curie temperature T c . 
P and F denote para- and ferromagnetic phases. Dotted line 
is a guide to the eye only. 

In summary, we presented a numerical proof of metallic 
ferromagnetism in the single band Hubbard model in the 
limit of infinite dimension. An asymmetric density of 
states, corresponding to a non-bipartite lattice, appears 
to be an essential ingredient for ferromagnetism. For a 
given interaction strength ferromagnetism is found over 
a wide hole doping regime but not for electron doping. 
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Finally, we investigate the doping dependence of the 
critical temperature, T c (Fig. 4). We find positive values 
of T c only in the case of hole doping and there over al- 
most the whole density regime. Only above a density of 
n 1=3 0.9 the extrapolation described above would lead to 
negative values of T c , but with the present approach we 
cannot rule out very low transition temperatures of the 
order of 10 -3 . For low densities the Curie temperature 
becomes very small and seems to vanish close to n = 0. 
Within the error-bars the results are consistent with the 
variational approach which yields a T c proportional to n 3 
for small n, i.e. no lower critical electron density [E4j. We 
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